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Abstract 



We develop a theory of existence, uniqueness and regularity for the following 
Oh ' porous medium equation with fractional diffusion, 

<' 

I ^ + (-A) 1 /2(| M p-i n ) =0 xeR N ,t>0, 
I at 

{ u(x,0) = f(x), x G R N , 

with m > = (N - 1)/JV, N > 1 and / G L 1 (R JV ). An L 1 -contraction 
semigroup is constructed and the continuous dependence on data and exponent 
is established. Nonnegative solutions are proved to be continuous and strictly 
00 \ positive for all x G R N , t > 0. 
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3 ■ 1 Introduction 

o ' 

This paper is concerned with the existence, uniqueness and properties of solutions 
u = u(x,t) to the Cauchy problem 

— + (-A) 1/2 (|u| m_1 n) =0, x e R N , t > 0, 
(1.1) < dt v ; Vl 1 ; 

{ u {x,0) = f(x), x e R N , 

for exponents m > 0, in space dimension N > 1, and with initial value / G L 1 (IR Ar ). 
By a solution it is meant a suitable concept of weak or strong solution. In particular, 
we prove that u G C([0, oo) : L 1 (IR Af )) and that the equation is satisfied a.e. in 
Q = M. N x (0, oo). The sign requirement u > is not strictly needed but when 
enforced some additional properties hold. 

We recall that the nonlocal operator (— A) 1 / 2 is denned for any function g in the 
Schwartz class through the Fourier transform, 

(1-2) (-~Ajv^(£) = i£l?(0, 
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or via the Riesz potential, 



:i-3) (-A)Wg(x) = C N P.V. [ t ] 7,ff dy, 



\x-y\ N+1 



JV+l . 



where — 7r 2 r( ; ^- i -) is a normalization constant, see for example 

Equations of this form can be considered as nonlinear variations of the linear frac- 
tional diffusion equation obtained for m — 1, which is a model of so-called anomalous 
diffusion, a much studied topic in physics, probability and finance, see for instance 
[H 123 [2H EH HH H2] and the references therein. We recall that fractional Laplacian 
operators of the form (— A)°"/ 2 , a E (0, 2), are infinitesimal generators of stable Levy 
processes [31 [11] . The analysis of the linear equation in the whole space is easy since an 
integral representation can be used for the solutions, see below. Such a representation 
is not available in the nonlinear case. 

Interest in studying the nonlinear model we propose is two-fold: on the one hand, 
experts in the mathematics of diffusion want to understand the combination of frac- 
tional operators with porous medium type propagation, and on the other hand models 
of this kind arise in statistical mechanics (27] and heat control [5] . The rigorous study 
of such nonlinear models has been delayed by mathematical difficulties in treating at 
the same time the nonlinearity and fractional diffusion. 

Observe that the above equation becomes the well-known Porous Medium Equation 
(PME) when replacing the nonlocal diffusion operator (—A) 1 / 2 by the classical Lapla- 
cian —A. A number of techniques in dealing with the present nonlinear fractional 
diffusion model will be borrowed from the experience obtained with the PME, see 
for instance [38J. Our original purpose was to study Problem (II. ip for every m > 1, 
to examine the existence and properties of "fractional slow diffusion". But the de- 
velopment of the theory allows to cover with a reasonable additional effort the "fast 
diffusion cases", m < 1, on the condition that we restrict the exponent to be larger 
than a critical value, m > m* = (N — 1)/N. This critical value is intrinsic to the 
equation, it appears in various contexts of the theory. It corresponds to the classical 
critical value m* = (N — 2) + /N in the PME case, see [9]. Existence of a weak solution 
is however proved for every m > 0. 

Harmonic extensions. Besides formulae fll.2p and fll.3p . there is another way of 
computing the half Laplacian, through the so-called Dirichlet to Neumann operator. 
If g = g(x) is a smooth bounded function defined in M. N , we consider its harmonic 
extension v = v(x, y) to the upper half-space v = EG?), i- e., the unique smooth 

bounded solution to 



(1.4) 
Then, 



v(x,0) = g(x), xeR N . 



;i.5) -^(x,0) = (-A x )^g(x) 

dy 
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where A x>y is the Laplacian in all (x, y) variables and A x acts only on the x variables 
(in the sequel we will drop the subscripts when no confusion arises). In order to 
check (11.51) . just apply the operator in the right-hand side twice. The choice of sign 
for the normal derivative makes the operator positive. Observe that the extension 
operator is well defined in H 1 / 2 ($L N ), and so is the Dirichlet to Neumann operator, 
which coincides with (— A) 1 / 2 in this more general setting. This well known technique 
has been recently used in several situations, see for instance [T6l [T7J [35] . 

Problem-setting. By means of the above-mentioned harmonic extension we 
rewrite, for smooth solutions, the nonlocal Problem (jl.ip in a "local way" (i.e., using 
local differential operators) as a quasi- stationary problem with a dynamical boundary 
condition. Indeed, w = |M| m_1 M satisfies 

{Aw = for x e R N , y > 0, t > 0, 

— - -A!_L- = for x G R N , y = 0, t > 0, 
w(x,0,0) = f m {x) for i G l w . 

This problem has been recently considered by Athanasopoulos and Caffarelli [5j . They 
prove that bounded weak energy solutions to (11.61) are Holder continuous if m > 1. 
The existence and uniqueness of that kind of solutions is one of the outcomes of the 
present paper. 

The connection between problems with dynamical boundary conditions and nonlocal 
equations has already been exploited in [10j in the case of a bounded domain, and 
in [2] for a semilinear problem in the half-space However, in those works the 

study of the nonlocal equation is used to obtain properties of the local one. Here, our 
approach is exactly the opposite. 

Main results. Our purpose is to establish a theory of existence, uniqueness, 
comparison and regularity for suitable weak solutions of Problem (11.11) with initial 
data / G L 1 {R N ). The full theory works for values of m larger than the critical value 
m* mentioned above, but basic existence and uniqueness holds for all m > 0, for data 
which are moreover bounded. 

Section [2] contains preliminaries, the basic definitions of solutions, and a list of main 
results. We define the concept of weak solution to Problem (ll.ip through the standard 
concept of weak solution to the associated local Problem (11.61) . We also define the 
concept of strong solution. 

We establish the existence of weak solutions in Section [3] by means of Semigroup 
Theory, solving first some associated elliptic problem, under the condition that the 
initial data / are both integrable and bounded. Actually, the obtained solution is 
strong and the equation is satisfied almost everywhere. We also prove in this context 
an L l -L°° estimate that will be basic for the so-called smoothing effect. 

Uniqueness is studied in Section [U Section \5\ deals with further properties of the 
constructed solution. It includes conservation of mass, positivity and regularity. At 
this point we use the continuity result from [5] to show that solutions to (11.11) corre- 
sponding to nonnegative initial data become immediately strictly positive if m > 1. 
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This is a remarkable property since it departs from the well-known properties of the 
standard PME, cf. [38]. On the other hand, for < m < 1 we are able to prove the 
expected positivity property using a different approach. This is used later, in com- 
bination with boundedness and a result in [5] to prove Holder continuity also in this 
case. Let us notice that, unlike in the local case, there is still no general regularity 
result for linear nonlocal equations (with reasonable coefficients) guaranteeing that 
positive bounded solutions to problem (11. ip are in fact C°°, though this is expected 
to be true. 

After such a work, we are able to treat general solutions with data in L 1 (M Ar ) in 
Section [6j Here we complete the proof of uniform boundedness of the solutions with 
integrable data for positive times, the L^-L 00 smoothing effect. 

In Section [7] we study the continuous dependence of the solution in terms of the 
exponent m and the data /, in the case m > m*. In particular, we show that the 
linear case m = 1 can be obtained as a limit of the nonlinear case both from above 
and below. 

Section [8] contains a brief description of alternative approaches to the existence 
theory and an announcement of extensions. Finally, the Appendix gathers some 
technical lemmas. 

Notice on the linear case. For the value of the parameter m = 1 we obtain the 
equation 



This is a linear fractional heat equation where the fractional derivatives act only on 
the space variable. It is explicitly solvable in terms of the initial value, u(x, 0) = f(x), 
through convolution with the explicit Poisson kernel in R^ +1 , 



where Cat is the constant in (II. 3ft . Note that this corresponds to an anomalous 
diffusion law of the form (x) ~ t a with a = 1 instead of the standard a = 1/2 of the 
Brownian case. 

Notations. In dealing with extended functions, we denote the upper half-space, 
by f2, and write its points as x = (x,y), x G M. N , y > 0. We denote by 
T the boundary of Q, i.e., F = M. N x {0}, which is identified to the original M. N 
with variable x. We consider also the extension and trace operators, E, Tr: for a 
function v G H l l 2 (T), we denote its harmonic extension to £1 as E(v); notice that 
Fi(v) G H l (Q); on the other hand, given a function z G if 1 (fi), we denote its trace 
on T, which belongs to i? 1 / 2 (r), as Tr(z). 

As in the PME theory, we will be mostly interested in nonnegative data and solu- 
tions. However, the basic theory can be developed for data of any sign, and in that 
case we will use the simplified notation u m instead of the "odd power" \u\ m ~ l u, and 



(1.7) 



du 
~dt 



+ (-A) 1/2 m = 0. 



(1.8) 
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we will also use such a notation when m is replaced by 1/m. In fact, we will show 
that if the initial value is nonnegative, then the weak solution we construct is also 
nonnegative, u > 0, which helps justify our abbreviated notations. 



2 Preliminaries and main results 

As mentioned above, we define the concept of weak solution to Problem fll.ip through 
the standard concept of weak solution to an associated local problem, which we write 
here again by convenience. 

Aw = for x G fi, t > 0, 

(2.1) { ^- — = onr.OO, 

w{x,0,0) = f m {x) onT. 

In order to define a weak solution of this problem we multiply formally the equation 
in (12.11) by a test function tp and integrate by parts to obtain 



T r rT 



(2.2) - / / (Vw, V(p) dxds + I I dxds = , 

'o Jn Jo 



with u = Tr(w 1/,,n ), on the condition that Lp vanishes for t — and t = T, and also 
for large \x\ and y. 

Definition 2.1 We say that a pair of functions (u, w) is a weak solution to Problem 
QUD if w £ V-^T^Wj^in)), u = Ti(w^ rn ) G L x (r x (0,T)) and equality (TO) 
holds for every tp G C^(Q x [0,T)). Finally, the initial data are taken in the sense 
that lim «(•,£) = f in L 1 (R N ). 

i — >0 

An alternative form of equality (I2.2p . including the initial value in it, is 



(2.3) 



f f - f f dt P 

/ / (Vw, Vip) dxds + / / u—dxds 

Jo Jn io ir ™ 



u(x, T)ip(x, 0, T) dx — I f(x)(f(x 1 0, 0) dx. 
I Jr 



As is usual, more general test functions can be considered by approximation, when- 
ever the integrals make sense. Note that the trace u = Tr(u> 1,/m ) is well defined. For 
brevity we will refer sometimes to the solution as only u, or even only w, when no 
confusion arises, since it is clear how to complete the pair from one of the compo- 
nents, u = Tr(w 1//m ), w = Fj(u m ). By a weak solution of our original Problem fll.ip 
we understand u, the first element of the solution to Problem (12. ip . 

Observe that the definite advantage of working with the local version is compensated 
in some sense by the difficulty of having integrals in (12. 2p defined in spaces of different 
dimensions. 
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This definition is a very general notion of solution: in this framework we can con- 
struct a weak solution to Problem (12. ip provided the initial value / is integrable and 
bounded. We restrict ourselves in the next results to such data. However, weak solu- 
tions are sometimes difficult to work with, and we are not able to prove uniqueness. 
Hence, a class of solutions with better properties is welcome. A quite convenient 
choice is the class of so-called weak energy solutions, cf. |38j for the standard PME. 

Definition 2.2 A weak solution pair (u,w) to Problem (II. ip is said to be a weak 
energy solution if moreover w G L 2 ([0,T]; H l (fl)). 

Theorem 2.1 Let m > 0. For every f G L 1 (R Ar ) n L°°(W N ) there exists a unique 
weak energy solution to Problem (12. ip . Moreover u G C([0, oo) : L 1 (M Ar )) nL°°(lR Ar x 



The importance of this class of solutions, besides having uniqueness, is that, if we 
restrict to nonnegative data and exponents m > m, = (N — 1)/N, we can obtain 
regularity and positivity. 

Theorem 2.2 Let f G L 1 (M. N ) n L°°(M. N ) be nonnegative, and assume m > m*. 
Then the weak energy solution (u, w) to Problem (12 .ip satisfies: 

i) Conservation of mass: for every t > we have 



ii) Positivity: u(-,t) > in M. N for every t > 0. 

Hi) Regularity: there exists some < a< 1 such that u G C a {R N x (0,T)). 

(iv) Maximum Principle: if U\,u 2 are solutions with data Uoi,«o2 and u m < u 02 a.e. 
in M. N , then ui < w 2 o- e., in Q = K x (0, oo). 

(v) Contraction: for any two solutions U\,Ui with data Woi,Wo2 we have 
(2.5) Wu^^t) - u 2 (-,t)\\ L i {R N ) < ||ttoi -«02||ii(RW). 

The restriction m > m* is not technical: positivity and conservation of mass are not 
true if m < m*, see Proposition 15.11 On the other hand, conservation of mass holds 
also for solutions with changing sign if m > m*. If m > 1 the C a regularity result is 
true also for any changing sign solution, [5]. 

A further interesting property is that the weak energy solutions are strong solutions, 
which means that the terms (in principle only distributions) involved in equation (12.11) 
are in fact functions, and equalities hold almost everywhere. The main technical 
difficulty is to prove that d t u is a function. 

Theorem 2.3 In the hypotheses of Theorem \2.2\ we have dtu G L 1 (M 7V ). 



[0,oo)). 



(2.4) 
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We observe that for strong solutions we can multiply (12. ip by any integrable function 
to get, instead of (12.21) . the following identity 



(2.6) / (Vw,V(p)dx+ / — (fdx = 0. 



Working with strong solutions we can use the solution itself as a test function in 
formula (12. 6p . In particular this allows us to obtain a universal bound for all solutions 
with the same mass. 

Theorem 2.4 Let f G L 1 (M. N ) nL oc (]R 7V ), and assume m > m*. Then, there exists 
a positive constant C such that the weak energy solution to Problem (12. ip satisfies 



(2.7) sup \u(x,t)\ < Cr 



j/N 



with 7 = (m — 1 + 1/AT) 1 . The constant C depends only on N and m. 

General integrable data. Once this theory is settled, we are interested in 
considering all integrable functions / as possible data in Problem (jl.ip . As we have 
advanced, this can be managed by approximation by bounded initial data, and this is 
possible if we have an /^-contraction at hand. We thus introduce the concept of L l 
energy solution: a weak solution, continuous in L 1 , which is also an energy solution 
for positive times. 

Definition 2.3 We say that a weak solution (u,w) to Problem (12. ip is an L x -energy 
solution if u G C([0, oo) : L 1 (M. N )) and \Vw\ G L 2 {Vt x [r, oo)) ; for every r > 0. 

The /^-contraction property for /^-energy solutions is as follows. 

Theorem 2.5 Let (u,w) and (u,w) be two L 1 -energy solutions to Problem (12. ip . 
Then, for every < t\ < ti, 

(2.8) / [u(x,t 2 ) — u(x,t2)]+dx < / [u{x, ti) — u(x, dx. 

Jr n ' Jr n 

We therefore have that, when performing the approximation by problems with 
bounded data, the limit function obtained is an L 1 -energy solution. Now, since es- 
timate (12. 7p does not depend on the L°° norm of the data, it is also true for the 
limit solution (for changing sign solutions it holds by comparison). In particular this 
represents an L 1 -L°° smoothing effect that allows to obtain the same properties of 
Theorem 12.21 for positive times. 

Theorem 2.6 Let m > m*. Then for every f G L 1 (M iV ) there exists a unique L 1 - 
energy solution to Problem ( 12. ip . It satisfies estimate (12. 7p and the conservation of 
mass (12. 4p . If moreover f > 0, positivity and regularity also hold, and the solution is 
strong. The maps S t : / u(t) generate a nonlinear semigroup of order-preserving 
contractions in L 1 (M. N ). 
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The next sections of the paper are devoted to treat the case of bounded data. In 
Section [6] we drop this restriction and deal by approximation with general L l data. 
We point out that the continuous dependence of the solutions constructed above with 
respect to the initial data and the exponent is stated and proved in Section [7J 

3 Weak solutions 

We set out to construct a weak solution to the extended local Problem ( 12.1 j) taking 
initial values / £ ^(T) R L°°(T). We point out that the construction of a weak 
solution can be made for every m > and data not necessarily signed. 

A well-known method of construction of solutions of evolution equations, and also of 
generating a semigroup in a convenient functional space, is the so-called Implicit Time 
Discretization. It runs as follows: if the evolution equation is dv/dt + A(v) = 0, where 
A is a linear or nonlinear, bounded or unbounded operator acting on a Banach space 
X, and given initial data v (0) = / £ X, then the construction of an approximate 
solution of the problem in a time interval [0,T] proceeds dividing the time interval 
[0, T] in n subintervals of length e = T/n and then defining the approximate solution 
v £ constant on each subinterval in the following way: in each interval (tk-i, = ke, 

k — 1, • ■ • , n, we consider the solution v e ^ to the discretized problem 



We take as starting condition v e q = / e , an approximation of /. In the case of 
linear operators, a variant of the Hille-Yosida Theorem ensures the convergence of 
these approximate solutions to the so-called mild solution of the evolution problem 
when the operator A satisfies some properties, like being maximal monotone, cf. [T3] . 
The convergence result in the case of nonlinear and possibly unbounded operators 
is given by the famous Crandall-Liggett Theorem [22] under the assumption that 
A must be accretive and satisfy a certain rank condition. (Reminder: a possibly 
nonlinear and unbounded operator A : D(A) C X — > X is called accretive if for every 
e > the map I + eA is one-to-one onto a subspace R e (A) C X and the inverse 
TZ(e, A) := (I + eA)~ 1 : R £ (A) — > X is a contraction in the A'-norm. The precise rank 
condition that we will use is R £ (A) D D(A) for every e > 0). 

One of the typical examples of such theory is the standard PME posed on the whole 
space or on a bounded domain with homogeneous boundary conditions. The early 
work due to Benilan and collaborators, [6], drew attention to this important results, 
as well as the application to more general nonlinear diffusion-convection models. 

We will apply such a strategy to our evolution Problem ( 12. lj) . The discretized 
problem is: 





S 



Given / G L 1 ^) D L°°(T) and e > 0, to find w e = {u e ,i, • • • , w e , n } by solving 
for k = 1, ■ ■ ■ ,n the problem 

Aw £: k = in £7, 

(3.2) <( dw B , k 



dy 



u £ ,k - u e u X on r 



with initial value w £j0 = / on T. In each such step, u £ ^-i = Trf^*^"^) is 
known and u e ^ and w e ^ = E('U e ,fc) are the unknowns. 

The second equation in f)3.2p can be written as 
(3.3) u £yk + sA{u £)k ) = u e ,k_i, 

where the operator A : -D(A) C ^(T) — >■ -^ 1 (r) is defined as 

'0E(v w 



(3.4) = - Tr 



<9y 



with domain 

(3.5) D(A) = {ve L\T) n L°°(r) : G L 1 ^), |M| LOO(r) < ||/|| LOO(r) } 

This operator is nothing but the half-laplacian of the power m, A{v) = (— A) 1//2 t> m . 

3.1 The elliptic problem 

Therefore, in order to perform the plan we need to establish the solvability and 
properties of the elliptic problem 

{Aw = in Q, 

— — h ur /m = g on 1 . 

for all g G -^(r) fl L°°(r). As we have said before, the power u> 1//m in the boundary 
condition means |iu| 1 ' m-1 M; if w takes on some negative values. We will also prove 
that if g > then w > 0. A weak solution to this problem is a function w G W l( £ (fi), 
such that Ty^w 1 /™ 1 ) G -^ 1 (r), verifying 



(3.7) 



I (Vw, + [ w^ m ip - [ gcp = 
Jn Jt Jr 



for any (p G Cq(O). We have to prove existence of the solution w and contractivity 
of the map g h- )■ Tr^ 1 /" 1 ) in the norm of -^ 1 (r), which plays the role of X in the 
definition of accretivity. 

To prove this we perform an approximation substituting the unbounded domain 
fl by an increasing sequence of bounded domains flft (half balls), imposing zero 
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Dirichlet condition on the part of the boundary of the domain which does not lie on 
the hyperplane y = 0. The approximate problems are 

Aw = o in n R = n n b r , 

(3.S) { — = w l l m -g onr fl = 9n fl n{y = 0}, 

iw = onE R = on R n {y > 0} , 

where B R = B R (0). The concept of weak solution for a given datum g G L 1 (r J j) fl 
L°°(T R ) is analogous to the one given above (13. 7p . after changing the domains of the 
integrals into the corresponding bounded domains. 

Theorem 3.1 For every g G ^(T) H L°°(r) there exists a unique weak solution 
w G H 1 ^) to Problem such that Ti(w 1/m ) G L\T) nL°°(T). Moreover, if w 

and w are the solutions corresponding to data g and g, then 

(3.9) J [w 1/m (x, 0) - w 1/m (x, 0)] + dx < J [g(x) - g{x)} + dx. 

This in turn implies that if g > in V then w > in Q. Moreover, ||ii||L°°(r) — 
\\g\\L°°(r)- 

Proof. Step 1. We first prove that there exists a weak solution w G H 1 ^^ to 
Problem ( 13. 81) . This is done by solving the following minimization problem: 

To find a function w G H 1 ^^ minimizing 

t/ \ 1 f In 12 m C \ m+1 f 

J(w) = - vw H / w \ m — / wo. 

This functional is coercive, since 

J(w) > Ci||w|||i (nji) - CsII^Hhm^), 

which follows by using the Poincare inequality, Cauchy-Schwartz and the trace embed- 
ding. Moreover, coercivity then provides a bound for HwIIh 1 ^)) though it depends 
on R. 

Step 2. We now establish contractivity of solutions to Problem (13 .8p in L 1 (r^j). 
Let w and w be the solutions corresponding to data g and g. We claim that 

(3.10) / [w 1/m (x, 0) - w 1/m (x, 0)] dx < [ [g{x) - g{x)} + dx. 

This inequality follows easily if we consider in the weak formulation the test function 
if = p(w — w), where p is any smooth monotone approximation of the sign function, 
< p(s) < 1, p'(s) > 0. We get 

p'(w -w)\V(w - w)\ 2 + (w 1/m - w 1/m ) p(w - w) -/ (g - g) p(w - w) = 0. 
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Passing to the limit, we obtain 



/ (w 1/m - w 1/m ) + dx < / (g - g) sign(w - w) dx < / (g — g) + dx. 
Jr R Jr R Jt r 

In particular, under the assumption g > we have w(-,0) > 0. Moreover, w(-,0) G 
L 1 (r^j) r\L°°(Tji). Finally, since the Poisson kernel in the half-ball is nonnegative, we 
also conclude that w > in fl R . 

Step 3. In order to pass to the limit R — > oo in the case of nonnegative data, we use 
a monotonicity property of the family of approximate solutions, denoted here by Wr. 
Namely, R < El implies wr < Wri in Qri. This follows from the ordering of the 
restrictions, using again that the Poisson kernel is nonnegative. The ordering of the 
restrictions results from comparison in r^, since Wr> > in E# (there is a contraction 
property analogous to f)3.10p for problems with non- homogeneous boundary data). 

Monotonicity implies that there exists the pointwise (and also in the sense of distribu- 
tions) limit w = lim wr. This limit satisfies w > in f2, w 1 / m (-, 0) G L 1 ^) nL°°(r). 

Since |Vwr| is uniformly bounded in L 2 (Qr), 



/ \Vw R \ 2 < / gw R < ||s'||li(m^)||5 , | 



we conclude that Vw^ — ^ Vw in L 2 (Q). This is enough to pass to the limit in the 
identity 

/ (Vw R , Vtp) + / w R /m tp - g^ = 0, 

2 



to show that w satisfies (13. 7L Also the estimate of the L norm of the gradients 
passes to the limit, and leads to 



(3.H) / |Vw| < \\g\\ L HRK)\\g 



Step 4. The pass to the limit in the case of non-positive data uses a similar argument. 
Finally, in the case of data g of both signs, we use comparison with the solutions with 
data gi = g + > and gi = —g~ < and compactness to pass to the limit. 

Step 5. Contractivity for the limit problem is proved exactly in the same way as 
for the approximate problems. This gives uniqueness. We also have that the L l norm 
and the L°° norm of the function w 1//m (-, 0) are bounded respectively by the L 1 norm 
and the L°° norm of the datum. □ 



3.2 Existence of solution for the evolution problem 

We now use the previously mentioned procedure to construct the solution to the 
evolution problem (12. lft . We recall that the Crandall-Liggett result only provides us 
in principle with an abstract type of solution called mild solution. 
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Theorem 3.2 For every f G -^ 1 (r) R L°°(r) there exists a weak solution (u,w) to 
Problem (EE]) withu(-,t) G L 1 (r)nL°°(r) /or every t > and w e L 2 ([0, T]; tf 1 ^)). 
Moreover, the following contractivity property holds: if (u,w), (u,w) are the con- 
structed weak solutions corresponding to initial data f, f , then 

(3.12) J[u(x,t)-u(x,t)] + dx < J[f{x) - f(x)] + dx. 

In particular a comparison principle for constructed solutions is obtained. 



Proof. For each time T > we divide the time interval [0, T] in n subintervals. 
Letting e = T/n, we construct the function w £ piecewise constant in each interval 
(tk-i, tk], where tk = ke, k = 1, • • • , n, as the solutions to the discretized Problems 
( 13. 2p . For convenience we write here again the problems: w £) k solves 

Aw £ ^ = in fl, 

dw £tk 

— = w £ ,fc - u e k-i on 1 , 

dy 

with u £ fl = f. Our solution is the (uniform in [0,T]) limit 

(u,w) = \im(u £ ,w £ ) 

in L} oc {Q). It is a mild solution, whose existence is guaranteed by the classical 
semigroup approach, u G C([0, oo) : L 1 (r)). In fact we obtain first the function 
u by Crandall-Ligget's Theorem, and the harmonic extension w of u m coincides with 
\im.w £ . By construction we have w G L°°(Q x [0,T]). We must now show that we 

have obtained in fact a weak solution. 

Multiplying the equation by w S) k, integrating by parts, and applying Young's in- 
equality, we obtain 



(3.13) e JJ\7w £ , k \ 2 dx< (m+1) ( / \u £ , k -i\ m+1 dx - / \u £tk \ m+1 dx) . 

Adding from k = 1 to k = n we get 



/ [ \VwJx,i)\ 2 dxdt < - — - — r / 
Jo Jn ~ (rn+l)J r 



\f{x)\ m+l dx. 



Passing to the limit, the same estimate is obtained for |Vu>|, and therefore w G 
L 2 ([0,T];/Z' 1 ( fi ))- 0n the other hand > fl3TT3|) yields 

\u £ A m+1 dx < J | M£ , fc _i| m+1 dx < J \f\ m+l dx, 

which, after passing to the limit, gives 

\u{x,t)\ m+1 dx< J \f{x)\ m+1 dx, 
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for every t G [0,T]. Now, choosing appropriate test functions, as in [31], it follows 
that we can pass to the limit in the elliptic weak formulation to get the identity of 
the parabolic weak formulation. We first have 

(Vw £tk , Vy?) = / -(w £|fe -i - u e>k )(p. 
Jr ^ 

Integrating in (tk-i, t k ) and adding for k — 1, • • • , n, we get that the right-hand side 
becomes 

n rt 



/ -(u £ ,k-i(x,t) - u £!k (x,t))Lp(x,0,t) dxdt 
fc=1 -Vi -/r ^ 



u £ (a;, t) -(^(a;, 0, t + e) — (p(x, 0, t)) dxdt 
o Jr £ 

+ - / / f(x)ip(x, 0, t) dxdt / / u £ (x, T)<f(x, 0, t) dxdt. 

6 Jo Jr e Jt-e Jr 

Passing to the limit e — > we get (12. 2p . 



The contractivity (I3.12p obtained in Theorem 13.11 in each step is inherited in the 
limit. In fact, if w £ ^ and w £ ^ are the discretized approximations of w and w, then 
we have 

J [u E)k +l{x) -U £)k+1 (x)} + dx < J [U e ,k(x) -Ue tk {x)) + dx, 

which easily implies (I3.12p . Comparison is a trivial consequence of contractivity. □ 

Remark. This contractivity property also implies the following estimates for the 
weak solution to Problem ( 11. ip just constructed 



(3-14) IK-,*)IUiCR*) - ll/IU 1 ^) ' \\u(-,t)\\ L ™ {U N) 



< 



Using now the Poisson kernel of the half-space and Young's inequality, we have that 
for every y > and every l/m<p<oo, it holds 

(3.15) HEK-x-,^)!!^) < o,t))|| LJW < H/lir^ll/ll^ 



with the first inequality replaced by equality if p = 1, m > 1. 

We end this section with a property satisfied by the weak solutions just constructed 
which is very useful in the proofs to come, with a number of other applications. 

Proposition 3.1 Assume f > 0. If w is the nonnegative weak solution to Problem 
(12. ip constructed in Theorem \3.2l then the inequality 

(3.16) (m-l)i^ + mffl>0 

at 

holds in the sense of distributions. 
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Proof. We use a simple scaling argument based on the homogeneity of the problem, 
(see jSj). Assume first m > 1. We have that for every A > 1 the function W\(x, t) = 
Xw(x, X^^t) is also a solution with initial value bigger than w(-,0). Then by the 
comparison principle we get, 

lim ^M)-^(M) 
a^i+ A — 1 

which gives (13.161) . For m < 1 the sign is reversed. Comparison can be easily justified 
in the discretized approximations. Since the function w\ is the limit of the rescaled 
approximations of the function w, we can compare w and w\. □ 

We observe for future reference that at the nonlocal level of function u we have the 
"monotonicity formulae" 

du u 

— > — - — it m > 1, 

(3.17) £* (m " 1}t 

— — < — it m < 1. 

<9t _ (1 - m)t 

Formula (13 . 1 6[) is empty for m — 1, but in this case it is easy to derive from the 
explicit representation of the solution that td t u + Nu > 0. 

Remark. In the PME model (the local analogue), a similar lower estimate of d t u is 
also available in the case m < 1. The proof uses in an essential way a second variable 
called the pressure, which is a potential for the velocity. It is not clear which could 
be the corresponding pressure for the nonlocal problem. 



4 Uniqueness of weak energy solutions 

In the previous section we have constructed a weak solution to the local Problem (12. ip . 
As we have said, the construction itself shows that this weak solution is in fact a 
weak energy solution. We prove uniqueness using an argument taken from Oleinik et 
al. [32]. 

Lemma 4.1 Assume f G L 1 (IR Ar ) fl L°°{R N ). There is at most one weak energy 
solution to Problem (II. ip . 

Proof. Let (u,w) and (u, w) be two weak solutions to Problem (12.11) . We take as 
test, in the weak formulation, the following function 

f T 

<p(x,t) = / (w — w)(x, s) ds, 0<t<T, 
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with ip = for t > T. Observe that this is a good test function when w and w are 
weak energy solutions. We have 



(V(iu — w)(x,t), / V(w — w)(x,s) ds) dxdt 



o Jn 

T 



(u - u)(x,t)(u m - u m )(x,t) dxdt = 0. 



'o Jr 

Integration of the first term gives 

•T 



1 
2 

+ 



V(w — w)(x, s) ds 







o Jr 



dx 



[u - u) (x, t) (u m - u m ) (x, t) dxdt = 0. 



Since both terms are nonnegative, they must be zero. Therefore, u = u in V. Obvi- 
ously this also gives, as a byproduct, w = w in Q. □ 



Remark. Observe that this proof only requires u(-,t), u(-,t) e L 1 (M Ar ) D L m+1 (R N ). 



5 Properties of weak energy solutions 

In the next sections we restrict ourselves to the range m > m*. 



5.1 Conservation of mass 

We establish next a property that is typical of diffusive processes. 

Theorem 5.1 If(u, w) is a weak energy solution to Problem (12. ip with initial datum 
f , then for every t > we have 

(5.1) / u(x,t)dx= / f(x)dx. 

Proof. The case m = 1 follows from the explicit representation (jl.8p . For general m we 
use the integral identity (12. 2p with a particular test function. Consider a nonnegative 
nonincreasing cut-off function ip(s) such that ip(s) = 1 for < s < 1, ip(s) = for 
s > 2, and define ¥>r(x) = i/j(\x\/R). We obtain, for every t 2 > h > 0, 

u(x, t 2 ) — u{x,t\) )(p R (x, 0) dx — — / / (Vw(x, t), Vipr(x)) dxdt. 
r v ' Jh Jn 

Integrating by parts, noting that — — (x, 0) = 0, we have that the space integral 

dy 

inside the right-hand side is 



1=1 w(x,t)A(fji(x)dx. 
15 



With the test function chosen we have 



\I\ < cR 2 w(x, t) dx. 

JR<\x\<2R 



>R<\x\<2R 

We now estimate this integral for every t\ < t < t 2 . If m > 1 we just observe that 

|/| <cR~ 2 / \w(x,y,t)\ dxdy < ci?" 1 t)||^" 1 RiV Ju(-, t)!^!^) 
Jo Jr n 
< cR' 1 -> , 

where we have used (I3.15p . In the case m < 1, applying Holder's inequality with 
some exponent p > 1/m we have 

|/| <cR- 2 \H-Ml~ ( 1 P N) ( I \w(x,t)\ 1 / m dx) 1/P \{R<\x\<2R}\^- 1 ^ 

V JR<\x\<2R ' 
-2R 



< C fl-^+ 1 )(cp-i)/rt||«(.,t)||^vp^ y^ y i^y.t)! 1 /"^)- 



Finally observe that if m > we can choose 1/m < p < N/ (N — 1) to force the last 
term to go to zero as R — > oo. □ 

Solutions that lose mass. We now present a result that shows the necessity of 
the condition m > m*. In fact if < m < m* there is a phenomenon of extinction 
in finite time, which makes impossible to have conservation of mass. The proof is 
almost exactly the same as the one in [H] for the PME model, where the corresponding 
condition on m is < m < (N — 2)/N instead of < m < = (N — 1)/N. 

Proposition 5.1 Let N > 1 and < m < (N - 1)/N, and let f G L X (M. N ) n 
L^~ m ' N (M. N ) . Then there is a finite time T > such that the solution u to Problem 
ffLTj) satisfies u(x, T) = in R N . 

Proof. As we have said, the proof follows the argument in [9]. Therefore we leave the 
details to be consulted there. Assume also for simplicity u > 0. 

If we consider = w p in the equality (I2.6P , we get 

4p f . e+1,3 , 1 [ d(u pm+1 ) , 

' \Vw—\ 2 dx+- / -dx = 0. 



{p+l) 2 Jn 1 + pm J T dt 

The use of this test function is justified in j9]. By the trace inequality we obtain that 
there exists a positive constant C = C(p, N) such that 

C^y u n-% dx) + J — — ^ — -dx < 0. 
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If we now choose p = (N(l — m) — l)/m > 0, we get {(p + l)mN) / (N — 1) = pm + 1. 
Therefore, the function J(t) = j T u pm+1 dx satisfies the inequality 

JSf—1 

J'(t) + CJ—(t) < 0. 

This implies extinction in finite time for J(t) and thus for u(x,t), provided J(0) is 
finite. □ 

Example. There exists an explicit example of the above extinction property for a 
particular m < m*. It has the form (separated variables), 

u(x,t) = G(x)H(t). 

Substituting this expression in (fTTTj) . we have H{t) = c(T - t) 1 ^ 1 '^, and G solves 
the nonlocal equation 

(-A) 1/2 G m = G. 

In the special case m = (N — 1)/(N + 1), there exists an explicit family of solutions 

G 6>T = A(r)[r 2 + \x-S\ 2 }- {N+m , 

with t > 0, 5 G and where A(r) > has an explicit expression, see [12J and 
also [19J. Observe that m < m* and G G L^IR^), and thus u(-,t) G L^M^) for any 
< t < T, while u(x,T) = 0. 

5.2 Positivity and regularity 

We show in this subsection that nonnegative bounded weak energy solutions are in 
fact positive everywhere in Q. This is true for every m > m*. The result is in sharp 
contrast with what happens for the local analog, the PME dtu = Au m in the case 
m > 1, for which initial values with compact support produce solutions that develop 
a free boundary. Free boundaries are a main feature of the standard PME theory, 
but they are not available here. 

The idea behind our positivity result is as follows: if u > is a classical solution and 
u(x ,t) = for some x and t, then formula (II. 3p gives (— A) 1//2 u m (x , t) < 0, and 
hence d t u(x ,t) > 0. If m > 1 we only have that the bounded solution u is Holder 
continuous by [5]. In the case m < 1 we do not even have that regularity. Now 
we perform rigourously the proof of positivity for weak solutions using the extension 
Problem (I21j) . 

Theorem 5.2 Assume m > m*. The weak energy solution u to Problem (II. ip with 
a bounded nonnegative initial datum is positive for positive times. Even more, the 
corresponding pair (u,w) to Problem (12. ip satisfies w(x,y,t) > for every x G M> N , 
y > and every t > 0. 
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Proof. The case m = 1 follows from formula (jl.8p . 

Case m > 1: We already know that for every t > 0, function w is Holder continuous, 
w > in fi, and u> > in Q. We have to prove that w > also on the boundary, T. 
By the comparison result, we only need to consider compactly supported initial data. 

In a first step we show that if w is not strictly positive on T at some time T, then 
the supports of w(-,0,t) form an expanding (in time) family of compact sets for 
< t < T. This follows from estimate (13. 17j) and Alexandrov's reflection principle. 
In fact (I3.17P implies that 

/tjNm/Cm-l) 
w(x,t 2 ) > W{X,ti)\^j- j , 

for every x E Q, t 2 > ti > 0, which is called retention property. Next we claim that if 
the support of the initial value / is contained in the ball Br, and w(xq, 0, t ) = for 
some xo G M. N , t > 0, then the support of w(-, 0, t) is also compact for any < t < t , 
and contained in a ball of radius depending on \x \ and R. 

To prove the claim we reflect, for any given point (xi,y) G fl, around the hyperplane 
7r = (xi — Xo) • (x — (xo + x\)/2) = 0. It is clear that if \xq — (xo + x\)/2\ > |x | + R, 
then the hyperplane ir divides the half-space Q in two parts, Q = Mi U M 2 with 
Br x [0, oo) C Mi, (xq, 0) G M-i. In this way, by the comparison principle, we obtain 
that the function z(x, y, t) = w(x, y, t) — w(xo + X\ — x, y, t) satisfies z(x, y, t) > for 
every (x, y) G Mi, t > 0. The comparison principle holds on the half-space Mi. Thus 
w(xi,0,t) < w(x ,0,t) = for every < t < t . A sufficient condition for x\ to get 
this argument work is |xi| > 3|xo| + 2R. 

In a second step we assume (thanks to the previous argument) that w(x, 0,t) = 
in some ball B C T for < t < t\. Then we have, for every test function ip that 
vanishes on dB x (0, oo), that 

f 11 f dip f h f°° f f h f dw 

= / / u ~^~ dxdt = / / (Viu, Vy?) dxdydt = — / / — — <p dxdt. 
Jo Jb ot Jo Jo Jb Jo Jb dy 

This gives d y w(x, 0, t) = for x G B, < t < t\. But w is a continuous nonnegative 
harmonic function in the half cylinder which vanishes on the part of the boundary 
y = 0. Hopf's Lemma implies d y w(x,0,t) < for x G B. This is a contradiction. 
Therefore, w(-,0,t) is positive everywhere. 

Case m < 1: In this case the proof is different, based on a weak Harnack inequality. 
Using estimate (13. 17ft and the fact that the solution is bounded, we know that for 
every t > there exists a constant A > such that 



(Vw, Vy?) + A wcp>0, 
n Jr 

for every nonnegative test function (p. Once we have this, we can use part of the 
proof of Lemma 2.4 in (15] to get a weak Harnack inequality in each large ball Br = 
{\x — xq\ 2 + y 2 < R 2 } C ~R N+1 with center on T. First of all, conservation of mass 



18 



plus (I3.15p . together with the fact that w > 0, w ^ 0, imply that there exists some 
R > large such that f D wdx > 0. 

° JtS R/2 

If we consider the function 

z(x,y) = e~ A ^w(x } \y\), 

then we have that z satisfies 

(Vz, V<^> -2A f sign(y)|^ -A 2 I z V > 0, 

i.e., it is a weak supersolution to an equation for which we can apply Theorem 8.18 
in [25] , to get 

inf z > cR-"- 1 1| z\\ L i {B 2) . 

£>R/2 

This means z > 0, (and thus w > 0) in Q. □ 

As a corollary of this result, we can establish also regularity for the case m < 1, 
provided m > m*. 

Theorem 5.3 Let m* < m < 1. T/ien any bounded weak energy nonnegative solu- 
tion u to Problem f 1 1.11) satisfies u G C" 1 ^ x (0, oo)) for some < a < 1. 



Proof. The above-mentioned regularity result of [S] applies to the equation 

+ (-A) 1/2 (t> ) = 



9 M , , aU/2, 



9f 

m some ball x <E B C R N , t > 0, with some nondegeneracy condition on the constitu- 
tive monotone function (3. Once we know that in such a ball the solution is bounded 
below away from zero, the requirements on the function (3 in [5] are fulfilled. □ 



5.3 Strong solutions 

We prove here that every nonnegative bounded weak energy solution is in fact a strong 
solution. We need to show that the time partial derivative of u is an L l function and 
that the second equality in (12.11) holds almost everywhere. 

As a first step we show that the time-increment quotients are bounded in -^ 1 (r), 
and thus the limit dtu must be a Radon measure. Our purpose is to prove that the 
limit is still in L l , and this is proved later. Observe that the result is clear in the case 
m — 1, from (jl.8p . We therefore assume 

Proposition 5.2 If u is the weak solution to Problem (11. ip constructed in Theo- 
remUJl then h~ 1 (u(-,t + h) - u(-,t)) G L 1 (R JV ) for every t, h > 0. 
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The proof is exactly the same as in the PME case, see [8]. In fact, the following 
estimate holds 



An analogous result for more general nonlinearities is given in [23] . 

Observe that by the Mean Value Theorem, the same type of estimate can be obtained 
for the time- increment quotients of w = u m if m > 1 since u G L°°. If m < 1 we 
obtain a bound locally for nonnegative solutions, since in each ball of ¥L N we have 
u > c> 0. 

We now prove a result which turns out to be fundamental in the proof of Theorem l2.6l 



(5.2) 



u(x, t + h) — u(x, t) 



Proposition 5.3 Let f > and let u be the weak solution to Problem ( II. ip con- 
structed in Theorem{2J\ Then <9 t (w (m+1)/2 ) G L 2 oc (R N x [0,T]). 

Proof. The formal proof is simple, using the function dtw as test function in the 
weak formulation, much as in the PME case, cf. [3_B1 Section 5.5], of course after 
some obvious changes. The problem is the justification of the calculations, since we 
have not established the existence of any kind of differentiability for the solutions or 
suitable approximations. 

Here we do by brute force as follows: 

• We use the weak formulation with the test function ip = 5 h (w* p), where p = rj*rj 
and r) is a convolution kernel, acting only on the time variable. We fix the notation 
f = f* r lif = f*Pi applied to functions of t. We also make use of the following 
calculus identity 

(5.3) J f(t)g(t)dt = j f{t)g(t)dt 

We prove this identity in Lemma [DTT1 of the Appendix for the reader's convenience. In 
addition, we take rj(t) = r] h (t) = (l/h)r]i(t/h), where rji is a smooth, symmetric, com- 
pactly supported nonnegative function, with support [—1,1]. Thus, p = ph inherits 
the same properties, with supp(pi) = [—2,2]. We also use the notation 

S h w(t) = 4" (w(t + h)- w(t - h) 
2h V 

for a discrete time derivative, omitting the spatial variable. Observe that since \Vw\ 
is in L 2 , this is a good test function. 

• Going to the weak formulation (12. 3p with the test function 8 w, we have 



(5.4) - f [ 5 h (w* p')udxdt+ [ 5 h (w*p)u\l 2 dx = - T f (V5 h w,Vw) dxdt, 
Jti Jr Jr 1 Jti Jn 

where p' = dp/dt. 
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The left-hand side mimicks ff(d t w)(d t u)dxdt = c JJ(d t (u^ m+1 ^ 2 )) 2 dxdt, while the 
right-hand side mimicks — jj (V(dw/ '<9t),Vii?) dxdt. Here the times t\ and t 2 are 
subject to be moved slightly on the condition that t\ < Ti < T 2 < t 2 for some fixed 
< Ti < T 2 . We now analyze the different integrals. 

• Using (15.31) . the integral in the right-hand side can be written as 
™*2 r 

(5 h Vw, Vmj) dxdt = 

'h Jn 

- 1 - f [ (Vw{t - h), Vw{t)) dxdt -^7 f [ (Vw(t + h), Vw{t)) dxdt = 

^ rti+h r -y rt 2 r 

— I / CVw(t - h),Vw(t)) dxdt / / CVw(t + h), Vw(t))dxdt = 

2hJ u Jn 2hJ t9 _ h J n 



2h Jt x Jn Jt 2 -h Jn 

— / / (Vw(t + h),Vw(t)) dxdt / / (Vw(t - h),Vw(t)) dxdt. 

'hJu-hJn 2h Jt, Jn 



2h 



1 ti —h J n Z ' IL J ti 

These two last integrals are bounded, in absolute value, by 

1 rti+h r 1 rt'2+h p 

— / / \Vw(t)\ 2 dxdt + — / / \Ww{t)\ 2 dxdt = Y l 
Jti-h Jn Jt 2 -h Jn 

Since |Vw| 2 G L 1 , by picking good times t± and t 2 from a diadic division of intervals 
around T\ and T 2 , the quantity Y\ is bounded uniformly in h (though ti and t 2 depend 
on h). 

• We turn to the left-hand side terms in (15 .4p . Assume first m > 1. The second 
integral in (15. 4p is bounded using Proposition 15.31 

As to the first integral in (15.41) we have 

t 2 r rtz—h /• 

/ 5 h (w * p')udxdt = / (w * p')5 h udxdt 

h Jr Jh+h Jt 

y r-U+h r- y rt 2 +h i- 

/ I ( w * P') u (t + h) dxdt — — \ (w * p')u(t — h) dxdt. 



2hJ tl _ h J r y r ' y ' 2h 



t 2 -h Jr 



The last two integrals are bounded in absolute value, in the same way as the second 
integral in (15.41) . The first term in the right-hand side above is the one we want to 
estimate carefully. By Lemma 19.21 and under the extra assumption of monotonicity 
in time, dtu > 0, we have 

(5.5) I = (w*p')(x,t) > c5 h w(x,t), 

with a positive constant c that depends on m, N and the smoothing kernel r]i. It only 
remains to observe that 

(5 h w) (5 h u) = (5 h u m ) (5 h u) > c(5 h u {m+1)/2 ) 2 . 

This turns out to be an easy calculus problem, using the technical Lemma 
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• Summing up, we get a uniform estimate 



f T2 I \5 h u^' 2 fdxdt<C. 
Jt± Jr 

Letting h 0, we get d^vr^ 1 ) G L 2 (IR Ar x [Ti,T 2 ]) for every < T x < T 2 . 

• General situation when m > 1: we want to apply Lemma [9.21 to I = (w * p')(x,t) 
without the extra assumption of monotonicity in time, only using that if m > 1, 
d t u > —cu. Then we use the second version of Lemma 19.21 An extra term appears 
but it is controllable. 

• In the case m < 1 we need to perform some little extra calculations. First, in order 
to use a bound of the form (15.21) for w we take advantage of the fact that w > c > 
in every compact set of Q. Thus we consider a new test function by multiplying ip by 
a cutoff function if)(x). The extra terms obtained are easily bounded. 

On the other hand, formula (15.51) holds with reverse inequality, 
(5.6) I = (w* p')(x,t) <c6 h w{x,t), 

provided d t w < 0, as before. Care has to be taken in the general case, where we use 
the estimate d t w < cw. □ 



We now prove the main result. 



Theorem 5.4 The weak (nonnegative) solution u to Problem ( II. ip constructed in 
Theorem 12.11 is a strong solution. Moreover, 



(5.7) 



du 



t) 



< 



\m 



It 



Proof. Let w be the weak solution to Problem (12. ip associated to u. We want to 
prove that the time derivative of u is actually an integrable function and that the 
normal derivative of w on T is a distribution. Thus the second equation in (12. ip . 
i.e. the equation (II. ip . holds almost everywhere. 

To deal with the normal derivative we only have to take into account that, thanks 
to the trace embedding, since for every t > we have w G H l (Q), then d y w(-, 0, t) is 
a distribution in H~ l l 2 (T). 

We now look at d t u(-, t). We use a technical result by Benilan f?j, see also [38]. As in 
the previous proof, we may assume dtu > 0. Also we know that u G L 1 (M Ar )nL 00 (lR Ar ) 
for every t > and finally, from Proposition ( 15. 3p we have dt(u^ m+1 ^ 2 ) G L 2 (R N ). 

All these estimates allow us to apply Lemma 8.2 in [38J to get d t u G Lf oc (R N ) for 
every p G [l,Pi), where p\ = min{(m + l)/m), 2}. This gives d t u G Lj oc (M, N ) for every 
t > 0. Formula (1577]) follows from (15/2]) . and in fact d t u G L 1 (R N ) for every t > 0. □ 

Once we know that our solution is a strong solution, we can perform the above- 
mentioned formal calculation, analogous to the case of the PME. 
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Proposition 5.4 In the above hypotheses, 
(5.8) 





Q u (m+l)/2 


n 


dt 



Tfi ~h 1 f 

(x,t)dxdt< / u m+1 (x,ti)cfe. 

Smti J r 



We end this section with two more estimates that will be useful in the sequel. By 
comparison assume / > 0. Putting (p = w as test function and obtain 

(5.9) [ [ \Vw\ 2 dxds + — '— / u m+1 dx = —^— [ f m+1 dx. 

Jo Jn m + lJr m + lJ r J 

We have thus a control of the L 2 norm of the gradient in terms of the initial data. 
Proposition 5.5 In the above hypotheses 

L m + 1 {vy 



(5.10) ||Viu|| J yi(nx(o I oo)) < c''^ 11 "" -1 - 1 '"' 



Another consequence of (15. 9ft is that the norm ||m(-, t) \\L m + 1 (R N ) is nonincreasing in 
time. In fact this also follows from the elliptic estimates of Section 13.21 

An easy generalization shows the following property. 

Proposition 5.6 In the above hypotheses any LP norm of the solution is nonin- 
creasing in time for every p > 1. Even more, if $ is any convex nonnegative real 
function, then J RN ^(u(x,t)) dx is a nonincreasing function. 

Proof. We have 

/ V(u(x,t 2 ))dx- [ V(u(x,t 1 ))dx = - I* - f w {1 - m)/m ^"(w 1/m )\Vw\ 2 dxds < 0. 
Jr Jr Jti m Jn 

□ 



6 Solutions with data in L 

In order to construct solutions with initial data in L 1 we approximate by problems 
with data in L 1 fl L°° and use the L 1 contractivity to pass to the limit (together with 
the estimates of the gradients in L 2 ). 

We begin by proving the bound in L°° in terms of the L 1 norm of the initial datum, 
estimate (12. 7\i . 
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6.1 Smoothing effect 



Proof of Theorem \2.4\ Case m > 1. We use a technique inspired by the proof in 
[26] . Lemma 3.3. To do that we consider the solution w to the local Problem (12. ip . 
Using estimate (13.171) and the integral formulation of that problem, identity (12. 6p . we 
get for any t > 0, 

/ (Vw, Vip) -L uip<0, 
Jn Jr 

where L = 1 /((m — l)t), <p > 0. The letter c will denote any constant depending only 
on m and N. Let (p = w p , p > 0. Manipulating the integral in Q and using the trace 
embedding, we obtain 



N-l 



Observe that 4p/((p + l) 2 ) < 1, so c does not depend on p. Therefore, if we denote 
\\u(-, t)||i?(r), for the fixed time t by ||w|| g , we have 



\\u\\ sm(p+1) < (cL)^\\u\\^ , 
where s = N/ (N — 1). We now iterate this estimate, 



\u\\ qk+1 <(cL)^\\u\\ q q k k + \ 



where q k +i = s{q k +m — 1). We obtain 



|u|| gh <(cL)°*|M|£, 



with the exponents 



1 h 



a k = — > s J , Ofc = 



If we start with q > 1, it is a calculus matter to obtain 

q k = As k -B, B = ^ -L = N(m-l), A = q + B, 

s — 1 

and then 

hm b k = — , hm a fc = — . 

We conclude that 

jv SQ. 
\\u\\oo < CLA \\ U \\ q A . 

Now recall that q > 1. In order to reach the L l norm, we use a classical interpolation 
argument once the above partial smoothing effect is obtained. To do this we put 
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go = 2 for simplicity, since any other value gives the same estimate. To prove an 
I}-L 2 smoothing effect we observe that 

IMh < \\u\\i \\ u \\o^ < cL2^a:||w||f \\u\\f*, 
where here A = 2 + N(m — 1). Therefore 

2 + iV(m-l) 



|| n || 2 < L 2(l+N( m -1)) 1^11^+^- 

which yields the estimate 



Recalling that L ~ 1/t, and using the Remark after Theorem 13.21 this gives (12. 7p . 

Case m < 1. Since we do not have a lower estimate for in this case, we use here 
a different idea: we look at the proof of the smoothing effect for the quasi-geostrophic 
equation performed in [18], see also [20]. In fact, the proof can be adapted to get also 
an L 1+£ -L°° re gularizing effect, for any e > 0. In order to reach L we use an analogous 
procedure as before, but here the interpolation must be performed iteratively, since 
the first regularizing effect in this case is obtained for different times. 

Let w be the solution to Problem (12. ip . and fix C > a constant to be chosen 
later. Consider the functions ( = (w — C) + and 9 = [w 1 /" 1 — C l / rn ) + . By Kato's 
inequality we have that ( is subharmonic. Using tp(0) as test function in the definition 
of Problem (12. ip , we get 



j t J i>(9) dx + J (V<P(8), VC> dx < 0, 



where ip' = ip. If we put <p{0) = 6 £ , the following will produce an L 1+£ -L°° regularizing 
effect. As we have said before, it is enough to get an L 2 -L°° regularizing effect in the 
first step. Thus we consider for simplicity f{6) = 9. We have then 

4- [ 2 dx + 2 [ (V0, V0 dx < 
at Jt Jn 

We now fix the constants. Let t > 0, t k = t (l - 2~ fc ), M > 0, C = C k = 
{M{1 - 2- k )) n \ 6 = 6 k ,C = Ck- Define the quantity, 

(6.2) U k = sup ( / 6 2 k {x, 0, t)dx] + 2 ! [ (V9 k , VCfc) dxdt = I k + Y k . 
t>t k \Jt J Jt k Jn 

Our purpose is to obtain a recursive estimate of U k in order to prove lim U k = 0. 

This will mean lim(Tr( M ; 1 / m )-C fc 1/m ) + = 0, i.e. u < M. We must choose M carefully. 

fe->-0 

First we deduce the upper estimate, for t k -i < s < t k , 

U k <2 J 9 2 k (x,0,s)dx. 
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Taking the mean over the interval [t k -x,t k ] we get 

2&+i r°° 



(6.3) 



o Jt k ^ Jr 



9 k (x, 0, s)dxds. 



See |l8j . To get a lower estimate, in the same way as in [IB], we have to take care of 
the exponents. Observe that the spatial integral in Y k can be written as 



8m 



I , , / m+1 

V(w 2m - C k 



m+1 
2m \ 12 



Using the trace embedding, and the fact that m < 1, we get 



m+1 



\7(w~ - C k 2m )+\ dx > c / \(w~ - C k 2m ) + | Ar - 1 dx 



m+1 



N-l 

"i+l 2iV ^ N 



1 J- iV(m + l) 

>c / |(^™-C fc m )+| ^-1 



jV-1 



iV(m+l) 

e k N ~ 1 dx 



N-l 
N 



where c depends only on m and N. From this estimate we obtain, using Riesz 
interpolation, 



u k > im- v > csu P ( / el 



> c 



t>t k 



N-l \ l—V 

JV(m + l) \ - j\T 



N-l 

iV(m + l) \ — — 



0^ 



1 








> c 


[f H 







l-v 



by taking z/ = 1/(N + 1), g = m + 1 + 2/N. Observe now that q > 2 provided 
m > 1 — 2/N, which is guaranteed by the condition m > m*. The last step is to use 
the above inequality to deduce the required iteration 



U k < 



2 k 



k+i r°° 



t. 



o Jt k _ t Jr 



0l(x, 0, s)dxds 



Ok+l+k(q-2) poo r r oKl-^) 

< - / / 0U(x, 0, s)dxds < C fj^Ul + J/ N . 



t Mi- 2 



t k -i 



This gives the limit mentioned above, lim U k = 0, provided 

k— >oo 



cU, 



1/N 



< 1. 



This can be achieved choosing M = ctQ^U^ N large enough, where l/^u = q — 2 = 
m — 1 + 2/N. This gives u < M. Since U < c||/||2, we obtain the L 2 -L°° estimate 



(6.4) 



\u(x, t)\ =< ct 
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We now perform the interpolation. Using ( 16. 4 p in the interval [t/2,t], and then in 
the interval [i/4, t/2], we get 

IK,f)lloo < c{t/2)-»\\u{.,t/2)\\f /N < c2^|| M (- ) t/2)||^ Ar || M (-,t/2)||^ 

<c2^|| M (- ) t/2)||^(c(t/4)^|| M (-,t/4)||^)" /7V . 

Iterating this calculation in intervals of the form [2~ k t, 2~( k ~ 1 't], using Proposition E21 
we obtain 

IK,0lloo<^2 6 *«-*||tt(. > 0)||?|K,2-*t)||?. 

Using the fact that m > m* implies /x < N, we see that the exponents satisfy, in the 
limit k — > oo, 

- (EV (rn-l)N + 2 _ 7 

(m-l)N+l~ N + ' 



fe— l 

$>(j+i)(^y -^6<cx- 

i=o 



dk = /iafc -> — — = 7 , 

[m — ljiv + 1 

7 

efc = a fc - 1 -> — , 

A = *(£)*^o. 

This ends the proof. □ 

Remark. Once we assert that the constant c in (12. 7p is universal, then the values of 
the exponents 7 and 7/iV are given as an immediate consequence of the invariance of 
the equation under the two-parameter scaling group. This is similar to what happens 
in the Sobolev inequalities, cf. [24J, page 262, or what happens in the PME, cf. [39] . 
page 29. 



6.2 L contraction 

Proof of Theorem \2.5[ We first recall that L^-energy solutions are weak energy solu- 
tions for every r > 0. Therefore, by the trace embedding we have u(-,t) G L P (M. N ) 
for 1 < p < 2N / (N — 1) if AT > 1 (for every p > 1 if N — 1). Now condition m > 
implies u(-,r) G L m+1 (R N ). This allows us to use the proof of uniqueness for weak 
energy solutions, cf. Lemma 14.11 and the Remark after it, for every r > 0. Conti- 
nuity in L 1 gives uniqueness up to r = 0. We now approximate the initial datum 
/ G L\R N ) by functions f k G L\R N ) H L°°(R N ), f k ->■ / in L\R N ). The above 
implies that corresponding solutions u k converge in L 1 (M 7V ) to our solution u, and 
the contractivity result in ( 13. 2 p passes to the limit. □ 
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6.3 End of proof of Theorem 12.61 



Once we have the L 1 contraction ( 12. 8p at our disposal, and using the above approx- 
imation of u by Uk, we observe that the estimate (12. 7p does not depend on the L°° 
norm of the approximations, and thus it is true for the limit. 

As to the conservation of mass, we recall that the proof of Theorem 15.11 relies on 
bounds for the L p norms of «(•, t) for t > 0. This is handled with the smoothing effect 
just proved. We thus have, repeating the proof of Theorem 15. 11 



u(x,t) - f(x))ip R (x) dx 



<cR~ 2 



\w(x, s) \ dxds. 



'0 J R<\x\<2R 

We now use (12.71) to estimate the right-hand side /. If m > 1 we get 

"/(m— 1) 



/ < cR- 1 



o 



u{-,s)\\ L i {R N)ds < cR 1 



i 



<r 7(m - 1} ds o, 



since 7(771 — 1) < 1. In the case m < 1 we obtain, instead, 

< c R N{p ' l)lp - 



7(m-l/p) 1 ft 
N ' p 



-l(m-l/p) ds ^ o 



by choosing 1/m < p < N/(N — 1), which is possible whenever m > m*. 

Finally, the proof of positivity and regularity in Theorems 15.21 and 15.31 also apply 
using again the smoothing effect. □ 



7 Continuous dependence 

The aim of this section is to prove the continuous dependence of the solutions con- 
structed in this paper with respect to the initial data / and exponent m. This is true 
for m> m*. Let us introduce the notation: if u = u m j is the solution corresponding 
to m and /, we write S(m, f) = u m j. 

Theorem 7.1 The map S : (m*, 00) x L\R N ) C([0, T] : L\R N )) is continuous. 



This will follow from a result of nonlinear Semigroup Theory which states that if each 
of A n , n = 1, 2, . . . , 00 is an m-accretive operator in a Banach space X, f n e D(A n ) 
and u n is the solution of 

— + A n u n = 0, u n {0) = f n , 
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then A n — > A^, f n — > implies u n — > in C([0, oo) : Af), where A n —> A means 
lim (/ + A n )~ x g = (I + A^^g for all g G X. 

n— >oo 

See, e.g, [21], [21] for statements and references. Hence, the theorem will be a corollary 
of the convergence of (/ + A^) -1 , where A m {u) = (— A) 1//2 w m = —d y w. This is what 
we prove next. 

Proposition 7.1 Lei {m n }^ =1 fre a bounded sequence of numbers in (m*,oo) such 
that lim m n = m>m t . T/ien lim (I + A mn )~ l g = (I + A fh )~ 1 g for all g G L 1 (R Ar ) fl 

n— >oo n— s-oo 

L°°(lR Ar ). 

Proof. We borrow ideas from [lOj. Note that an analogous result for the PME was 
proved in [9]. Let u m = (I + A m ) -1 ^. The L 1 -contraction estimate (13. 9 p implies the 
bounds 

ll M m n |Ul(RW) < IblU 1 ^) j 

(7.1) 

\\u mn - r h u mn \\ L imN) < \\g - r h g\\ L i 



for each h G M^, where (r/ l f)(x) = v(x + /i). This is enough, thanks to Frechet- 
Kolmogorov's compactness criterium [T3], to prove that {u mn } is precompact in L l (K) 
for each compact set K C R . 

To extend compactness to the whole we need to control the tails of the solutions 
uniformly. More precisely, we need to prove that, given e > 0, there exists some 
R > such that \\u mn \\ L i(^N\ BR ^ < e. This follows from a computation which is 
very similar to the one in the proof of Theorem I5.1[ but now taking as test function 
1 — <p R instead of (p R . First observe that u mn — (I + A mr ^)~ x g means 

/ u mn Lp R dx= I g<p R dx- / (Vw mnJ Vip R ) dx. 
Jr Jr Jn 

Therefore, with the above mentioned test function we obtain 



u mn dx < gdx + cR / w mn dx — > 

'{\x\>R} J{\x\>R} J{R<\x\<2R} 

as R — y oo uniformly in n, see the proof of Theorem 15.11 

We have obtained that along some subsequence, which we also call {m n }, the fol- 
lowing convergence holds 

u m „ ->■«*, in ^(r), 

for some function u*. Using the Poisson kernel, we also have w mn — > w* in L 1 ^), 
where w* is the harmonic extension of w* to the upper half-space. On the other hand, 
we have a uniform control in L 2 (Q) of the gradients of w in terms of the data g, 
see (13. lip . Thus, there is weak convergence in L 2 (Q) of the gradients Vw m „ along 
some subsequence towards Vw». All this is enough to pass to the limit in (17. ip to 
show that the limit is indeed u m . □ 



Remark. If m > 1, an easier alternative proof can be performed using the compact- 
ness results of |36|. 
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8 Comments and extensions 

Alternative approaches. It is not difficult to prove a posteriori that the con- 
structed semigroup is also contractive with respect to the norm H~ 1 / 2 . This property 
could be used as a starting point in the existence and uniqueness theory by using 
results of the theory of monotone operators in Hilbert spaces, as developed in [13] 
for the PME case. We have chosen our present formulation because we have found a 
number of advantages in proceeding in this manner. 



On the other hand, Crandall and Pierre developed in [23] an abstract approach to 

study evolution equations of the form d t u + A(p(u) = when A is an m-accretive 
operator in L 1 and (p is a monotone increasing real function. This allows to obtain 
a mild solution using the Crandall-Liggett Theorem. Our problem falls within this 
framework. Let us point out that such an abstract construction does not give enough 
information to prove that the mild solution is in fact a weak solution, in other words, 
to identify the solutions in a differential sense. 

Extension. As a natural extension of this work we can consider the more general 
model based on the equation 



where the fractional Laplacian has an exponent a G (0,2), and m > m* for some 
m*(N,o~) G [0,1). Though the main qualitative results are similar to the ones pre- 
sented here, the theory of these fractional operators with a ^ 1 has some technical 
difficulties that make it convenient to be treated at a second stage. We recall that 
Caffarelli and Silvestre [17] have recently characterized the Laplacian of order a, 
(—A) "/ 2 , by means of another auxiliary extension approach. We will use such an 
extension in a separate paper, [33], to treat in detail the more general fractional 
diffusion model. 

9 Appendix 

We recall some technical results that we have needed in the proof of the property of 
strong solutions. The first seems to be well-known. 

Lemma 9.1 On the condition that 77 is a smooth convolution kernel with r](—x) = 
rj(x) we have for every pair of L 2 functions in R 



Lemma 9.2 i) Let g be a positive nondecreasing function. If Ph(t) = (1 / h) pi(t / h) , 
where p\ is a smooth, symmetric, compactly supported nonnegative function, with 




(8.2) 



du 
dt 



+ (-A) CT / 2 (| M | m - 1 M ) = 0, 



(9.1) 
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support [-2, 2], and <5^(t) = (g(t + h) - g(t - h))/(2h), then 

(9.2) (g*p' h )(t)>c5 h g(t). 

ii) If instead of nondecreasing we have g'(t) > —Ag{t), the conclusion is 

(9.3) (g*p' h )(t)>c5 h g(t)-cAg(t). 



Proof, i) In the case of nondecrasing g we do as follows: 

(g*p' h )(t) = J (g{t + s)-g{t-s)){-p'{s))ds 



> 



f K (g{t + h)-g{t-h)){-p'{s))ds 

>5 h g(t) / (-hp'(s))ds = 5 h g(t) Pl (l). 
Jh 

We pass from the fist to the second line using the positivity of the integrand and the 
fact that for s e (h, 2h) we have g(t + s) > g(t + h) > g(t — h) > g(t — s). 

ii) The difference is now that we have to use the weaker inequality g(t + s) > g(t)e~ As 
if s > 0. We have again 

(g * p'hW) = J o (g(t + s)- g(t - s)) {-p'(s)) ds. 

Since g(t + s) - g(t - s) > -g(t)(e As - e~ As ) we get 

£ (git + s) - g{t - s)) (-p'(s)) ds > -g{t) f\e As - e- As )(-p'(s)) ds 

and the last integral is bounded uniformly for h small in the form cA. For the other 
part of the integral we have 

J {g{t + s)-g{t-s)Y-p\s))ds 

> / (g(t + h)e- A ^ -g(t-h)e A ^j(-p'(s))ds 

> 5 h g(t) / (-hp'(s)) ds-h = S h g(t) Pl (l) - h , 

Jh 

where h = (g{t + h)(l- e~ A{ - s -^) + g(t - h)(e A( < s - h *> - 1)) (-p'(s)) ds can also be 
estimated as cAgit). □ 

We end this list of results with an easy but useful calculus lemma. 
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Lemma 9.3 There exists a positive constant c depending on m > such that 

(9.4) (x m -l)(x- 1) >c(x^ -l) 2 , Vx>l. 
and also 

(9.5) (x m + l)(x + l) >c(x^ + l) 2 , Vx>l. 



Proof. The quotient of the two positive functions F(x) = (x m — l)(x — 1) and 
G(x) = (x^2~ — l) 2 is bounded below away from zero in the interval [1, oo) unless the 
limit at x = 1 or x — > oo is zero. It is clear that at infinity the limit is 1, whereas at 
x = 1 we can use L'Hopital to get 

F(x) Am 
i™ G{x) = (m + l) 2 ' 

This number is positive and equal or less than 1. The other inequality is similar. □ 
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